addpath('npy-matlab-master\npy-matlab')
A1 = readNPY('map_MCD_30K_niI2_test_y.npy')
surf(A1,'edgecolor','none')
set(gca,'xlim',[1 41],'ylim',[1 41])
view(2)
colormap(winter)
caxis([2.5, 6])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
%colorbar()
set(gcf,'Position',[0, 0, 300, 300])